Eliashberg theory of superconductivity and inelastic rare-earth impurity scattering in 

filled skutterudite La 1 _ 2 .Pr :I .Os4Sb 12 



Jun Chang 1 , Ilya Eremin 1 ' 2 , Peter Thalmeier 3 , and Peter Fuldc 1 
1 Max-Planck-Institut fur Physik komplexer Systeme, D-01187 Dresden, Germany 
2 Institute fur Mathematische und Theoretische Physik, 
TV '-Braunschweig, D-38106 Braunschweig, Germany 
3 Max-Planck-Institut fur Chemische Physik fester Stoffe, D-01187 Dresden, Germany 

(Dated: February 2, 2008) 

We study the influence of inelastic rare-earth impurity scattering on electron-phonon-mediated 
superconductivity and mass renormalization in Lai-^Pr^C^Sb^ compounds. Solving the strong- 
coupling Eliashberg equations we find that the dominant quadrupolar component of the inelastic 
scattering on Pr impurities yields an enhancement of the superconducting transition temperature T c 
in LaOs4Sbi2 and increases monotonically as a function of Pr concentration. The calculated results 
are in good agreement with the experimentally observed T c (x) dependence. Our analysis suggests 
that phonons and quadrupolar excitations cause an attractive electron interaction which results in 
the formation of Cooper pairs and singlet superconductivity in PrOs4Sbi 2 . 

PACS numbers: 74.70.Tx, 74.20.-z, 71.27. +a 



The Pr-based filled-skuttcruditc compounds have at- 
tracted much attention because of their exotic prop- 
erties like metal-insulator transition or unusual heavy- 
fermion behavior^i^ This also concerns recently dis- 
covered PrOs4Sbi2, the first Pr-based heavy- fermion su- 
perconductor with T c = 1.85 K which possesses many 
exotic properties compared to that of the Ce- and U- 
based superconductors 4 The heavy-electron mass has 
been confirmed by the large specific heat jump AC/T C ~ 
500 mJ/ (K 2 mol) at T c (see Ref. 0) and by de Haas- van 
Alphen (dHvA) measurements^ The ground state in the 
crystalline electric field (CEF) for Pr 3+ ion is a Ti sin- 

glet, which is separated by the first excited state, a T\ 
triplet by a gap of Acef ~ 8 K£ Because of the small 
Acef-, the relation between the quadrupole fluctuations 

(2) 

associated with the T\ state and the superconductivity 
has been recently the focus of intense discussions^ 

At present, the Cooper-pairing mechanism and the 
corresponding symmetry of the superconducting gap in 
PrOs4Sbi2 is still under debate. For example, initial 
studies of the thermal conductivity in a rotated mag- 
netic field have suggested the presence of two distinct SC 
phases^ Similarly, the London penetration depths has 
indicated a possible nodal structure of the superconduct- 
ing gap. On the contrary, a number of experimental tech- 
niques such as nuclear-quadrupole-resonance (NQR) ; 11 : 12 
Scanning Tunneling Microscopy (STM)^ muon spin re- 
laxation, (fj,SR)H thermal conductivity^ and specific 
heat measurements^ are in agreement with a fully de- 
veloped isotropic s-wave-likc superconducting gap at 
the Fermi surface (FS). Another interesting experimen- 
tal observation by means of the dHvA effect was that 
PrOs4Sbi2 has a very similar FS topology as the conven- 
tional s-wave superconductor LaOs4Sbi2 (Ref. [H) (T c = 
0.74 K) which hints at a relatively weak hybridization 
of the conduction band with Pr 3+ 4/ 2 electrons. Fur- 
thermore, it has been foun d 12 ' 17 that in the alloy com- 
pound Lai_ 2; Pr T Os4Sbi2 the superconducting transition 



temperature changes smoothly upon changing x. These 
observations raise doubts about the unconventional na- 
ture of the Cooper pairing in PrOs4Sbi2- But the en- 
hancement of T c and the effective quasiparticle mass in 
this compound with respect to that of LaOs4Sbi2 have 
to be understood. 

It is well known that scattering by magnetic impurities 
suppresses conventional s-wave superconductivity by de- 
stroying Cooper pairs in a singlet state. This situation, 
however, can change in the case of paramagnetic non- 
Kramers rare-earth impurities. For example, it has been 
shown previously^ that for superconductors containing 
an impurity with crystal-field-split energy levels the in- 
elastic charge scattering of conduction electrons from the 
aspherical part of the 4/ charge distribution may yield 
an increase of T c . However, in most of the cases the 
usual magnetic exchange interaction is dominant, thus 
suppressing T c . In this Rapid Communication, we solve 
the nonlinear Eliashberg equations for electron-phonon 
mediated Cooper pairing including inelastic scattering 
on Pr impurities. We find that a dominant quadrupo- 
lar scattering introduced by the Pr ions is responsible for 
an increase of T c in Lai_ x Pr2;Os4Sbi2 as a function of Pr 
concentration. 

Before solving the nonlinear version of the Eliashberg 
equations let us illustrate the results of Ref. [ID where 
the influence of inelastic scattering by rare-earth impu- 
rities on the superconducting transition temperature T c q 
without impurities has been studied for both exchange 
and quadrupolar scattering. In particular, for impurities 
with two singlet levels separated by an energy Acef the 
analog of the Abrikosov-Gor'kov relation for the change 
of T c as a function of the total impurity scattering ratei^ 
for the s-wave superconductor is given by 
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FIG. 1: Calculated values of the superconducting transi- 
tion temperatures as a function of Acef /T c q for inelastic 
quadrupolar (a) and magnetic (b) scattering by impurities 
with two singlet levels using Eqs. (1). 
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tcring between T 4 and T\ will be quadrupolar as well. 
Furthermore, from Fig.QJa) we see that the enhancement 
of T c is largest for Acef ~ 8 K. Keeping in mind that the 
higher excited states lie at energies > 100 K the magnetic 
pair-breaking scattering from these levels will play only 
a minor role in affecting T c . This suggests that Pr 3+ ions 
in Lai_ x Pr I Os4Sbi2 yield a strong Cooper-pair enhance- 
ment as compared with LaOs4Sbi2. At the same time, 
the Pr 3+ ions in Lai_ a; Pr. I ,Os4Sbi2 can be treated as inde- 
pendent impurities only at low enough doping. At large 
Pr concentrations the RKKY interaction between the 4/ 2 
ions results in the formation of O^y-type quadrupolar 
exciton bands as confirmed experimentally^ The latter 
can be treated as low-energy bosonic excitations which 
contribute in addition to phonons to the Cooper-pairing. 
But one should notice that the exciton dispersion is rel- 
atively weak, i.e., about 10% of Acef~ 

Using the generalized Holstein-Primakoff (HP) 
method^ the low-energy Hamiltonian describing the 
magnetic and quadrupolar interaction between conduc- 
tion electrons and the CEF split energy levels of Pr 3+ is 
given by 
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Here, 1/rfi and 1/r-g are the magnetic and quadrupo- 
lar scattering rate, respectively, and x = Acef /^ksTco. 
Furthermore, A(x) and B(x) are the combinations of 
digamma functions 1 - and arc smooth functions of x. 
In Figs. Ola) and IHb) we show the behavior of 
T c /T c q as a function of Acef I^bT c q for the quadrupo- 
lar (charge) and exchange (magnetic) scattering, respec- 
tively. Assuming equal matrix elements for both scat- 
terings (rf^Tco = t® 2 T c o = Ti2T c o) and setting some- 
what arbitrarily tyiT c q = 1.4 we find very different be- 
havior. While for quadrupolar scattering T c increases 
with respect to T c o, with a pronounced maximum around 
AcEF/ksTco ~ 10, the magnetic scattering results in a 
complete suppression of T c at small Acef and gives a 
moderate reduction at larger x values. Overall one finds 
that for similar scattering amplitudes the magnetic ex- 
change scattering dominates and T c decreases. 

Let us consider more specifically Pr 3+ impurities in 
Lai_ 2; Pr a ;Os4Sbi2 compounds. As mentioned above the 
CEF level scheme is composed of a Y\ singlet ground 

(2) 

state, a T\ ' triplet first excited state with Acef ~ 8 
K, and other excited states located at much larger en- 
ergies above ~ 100 K.— Furthermore, due to Th crystal 

(2) (2) 

group symmetry, the T y 4 ' level is of the form \T\ (m)) = 
VI - d 2 \T 5 {m))+d\T 4 (m)), where m = +,-,0. Here, T 5 
and T4 are the wave functions for the cubic Oh symmetry 
and d is a parameter that describes the mixture of the two 
states and is estimated to be equal 0.26J£ Since scatter- 
ing between T\ and Ts is quadrupolar, the dominant seat- 



where w q is the energy dispersion of the exciton. There 
are three bosonic modes (u = a, b, c) corresponding to the 

(2) 

excitations between singlet Fi and triplet T\ states. <x 
are the Pauli matrices (c = x, b = z, a = y) and / c (q) = 
q x q y , /b(q) = and / Q (q) = q y q z arc the form fac- 

tors of the quadrupolar interaction. The quadrupolar 
and magnetic excitations are written in terms of the 

bosonic operators O qu = i f/3 qu — /3l qM J and J q[t = 
(/3qu + /3l qu ) with (A«) 2 -(A cef + 2D M zjq)/ojq and 

(A q 7 ) 2 = (Ac7_E_F + 2L»QZ7 q )/w q , respectively. Here, D M 
and Dq are the effective magnetic and quadrupolar cou- 
pling constants^ z is the coordination number, and 7k 
is the corresponding structure factor. Due to the weak 
dispersion of the excitons, A^j and A q f arc almost momen- 
tum independent. Solving the Eliashbcrg equations, we 
use an effective interaction between quadrupolar excitons 
and conduction electrons averaged over the FS. 

On the real frequency axis the finite tempera- 
ture Eliashberg equations for the superconducting gap 
A(w,T) and the renormalization function Z(u>,T) are 
given by2i 
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and #±(fi) = a£G(n) + c^F(fi) =p a 2 M F(fl) is the gen- 
eralized electron-boson coupling function averaged over 
the FS. Here, n* is the screened Coulomb repulsion, f(ia) 
and n{ui) are the Fermi and Bose distribution functions, 
respectively, and /3 = l/feeT. 

As mentioned above there are two types of bosonic ex- 
citations. In LaOs4Sbi2 the superconductivity is driven 
by the electron-phonon interaction and we assume a sin- 
gle Lorentzian phonon mode at Tilue = 26 meV with 
broadening T p = loe/5. This corresponds to the De- 
bye temperature of this compound^ (see Fig. [5]). Setting 
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FIG. 2: Calculated bosonic Eliashberg function H±(Q) for 
Lai_ a ,Pr a; Os4Sbi2. The solid curve denotes the phonon con- 
tribution while the dashed and dotted curves refer to the 
quadrupolar and magnetic contributions, respectively. 



IX* = 0.1 and A = 2 f_ dna 2 P G(fl)/n = 0.33 the Eliash- 
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berg equations yield T c = 0.74 K and 2A /fcsT c « 3.5. 
In Fig. [3{a) we show the calculated results for the real 
part of A(u>) and the renormalization function Z(lo). In 
both cases one finds typical behavior of a BCS-like super- 
conductor. In particular, ReA(u;) shows a peak structure 
at uje and becomes negative at larger frequencies, reflect- 
ing the effective repulsion for lo > cue- For loe >> Aq, 
the renormalization function is close to 1 and the ratio 
2Ao/A:bT c equals the BCS value. 

A similar crystallographic structure and nearly equal 
ionic radii of the La and Pr ions allow to assume nearly 
the same phonon modes in LaOs4Sbi 2 and PrOs4Sbi2. 
The main difference is the occurrence of the low-energy 
exciton in PrOs4Sbi2 at about Acef ~ O.OSloe- As dis- 
cussed above there are both (magnetic) pair-breaking and 
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FIG. 3: Calculated frequency dependence of the real part 
of the superconducting gap function A(lo) for LaOs4Sbi2 (a) 
and PrOs4Sbi2 (b) compounds. The renormalization func- 
tion, Z(u>) at T = 0.04 K is shown in the insets. Note that we 
set the cutoff frequency equal to 5loe and introduce a finite 
damping T = 0.01 meV. 



(quadrupolar) pair-forming scattering processes of elec- 
trons by the excitons. Both processes can be taken into 
account as separate contributions to the Eliashberg func- 
tion a 2 F(Lo). Using Eq. §2fy and ignoring a weak disper- 
sion of the exciton as observed by inelastic neutron scat- 
tering (INS) experiments^ we obtain, after averaging over 
the FS, A Q = 2 r °° dLoa 2 Q F(co)/Lo a I 2 ac \M^\ 2 2N / A CEF 
where |M^| 2 = 35(1 — d 2 ) is the quadrupolar scatter- 
ing matrix element and N = 7.625 eV _1 is the bare 
(unrenormalized) electronic density of states^ Further- 
more, we set I ac « 2 meV. Similarly, we find A a/ = 
2j oo du J a 2 M F(L0)/u J » (1 - g L fll\M*i\ 2 2N /A C E F 
where, M$ = 20d 2 /3 is the magnetic scattering ma- 
trix element and gL =0.8 is the Lande' factor. Following 
our previous estimates, we set ag/a^ = 7 and obtain 
I ex sa 16.4J oc , where I ex is the on-site magnetic interac- 
tion constant between exciton and conduction electrons. 
All components entering the total Eliashberg function 
H±(Q) are shown in Fig. [2] As expected the quadrupo- 
lar contribution is about 7 times larger than the mag- 
netic one. Therefore an inclusion of the excitons yields 
an enhancement of the superconductivity in PrOs4Sbi2 
as compared to LaOs4Sbi2. In particular, the total value 
of A increases to 3.05 yielding strong coupling supercon- 
ductivity in PrOs4Sbi2 with T c = 1.85 K. Here, one has 
to remember that in PrOs4Sbi2 there are three contri- 
butions to the total A and their sum gives A=3.05. This 
yields a strong renormalization of the quasiparticle mass 
in PrOs4Sbi2 compared to its LaOs4Sbi2 counterpart. 
At the same time T c and A(lo) are determined by Eq. (5) 
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where the difference between quadrupolar and magnetic 
scattering enters. Therefore the enhancement of T c is 
moderate despite the large change in A. 

Figure[3]Jb) shows the calculated frequency dependence 
of the real part of the superconducting gap and the 
renormalization function for PrOs4Sbi2- One notices a 
sharp feature at an energy of about Acef + Aq and a 
strong renormalization of the quasiparticle mass due to 
the presence of the low-energy exciton. It is interest- 
ing to note that both magnetic and quadrupolar scat- 
terings contribute to the renormalization of the quasi- 
particle mass while only their difference contributes to 
Cooper-pair formation. We suggest that the low-energy 
peak in A(lo) should be visible in the tunneling density 
of states which could be checked experimentally. An- 
other remarkable feature is that we find 2Ao/fcsT c «5.4 
typical for a strong-coupling superconductor. Therefore, 
within this approach we find PrOs4Sbi2 to be a con- 
ventional s-wave strong-coupling superconductor, with 
enhanced T c due to low-energy excitons with mainly 
quadrupolar scattering. It has been proposed previously 
that quadrupolar scattering may alone yield unconven- 
tional superconductivity in PrOs4Sbi2 due to the exci- 
ton dispersioni^ Indeed such a possibility cannot be ex- 
cluded from our calculation. At the same time we find it 
difficult to reconcile with the experimental evolution of 
T c in Lai_ x Pra;Os4Sbi2 as a function x. In particular, 
assuming the linear dependence of the exciton's contri- 
bution to H±{lo) as a function of Pr concentration we 
find good agreement between calculated and measured 
T c (x) [see Fig. Ufa)]. At the same time, unconventional 
superconductivity in PrOs4Sbi2 would be unstable with 
respect to nonmagnetic La impurities and T c should van- 
ish at a small concentration of La. In Fig. Q] (b) we show 
the evolution of the effective mass as a function of Pr 
concentration. We find that the effective mass increases 
by a factor of 2.5 from LaOs4Sbi2 to PrOs4Sbi2 which 
agrees with dHvA data£ An inclusion of the higher-lying 
CEF energy levels provides stronger renormalization of 
the effective mass.— It is interesting to note that ex- 
perimentally AC'/T C shows a nonlinear dependence on 
the La concentration, in particular close to PrOs4Sbi2 17 
which cannot be explained by our theory. Whether this 
is a sign of the first-order phase transition between differ- 
ent symmetries of the superconducting order parameter 
or is a result of the Pr-Pr interaction neglected in our 
theory remains to be understood both theoretically and 



experimentally. Finally, let us note that our theory al- 
lows us also to explain qualitatively the behavior of T c 
in Pr(Osi_ x Ru x )4Sbi2 for x < 0.6. In particular, it has 
been found that the CEF level splitting between the low- 
est levels increases as a function of x while T c decreases^ 
According to our Fig. []Ja) the Cooper-pair constructive 
quadrupolar scattering decreases with increasing x while 
the magnetic pair-breaking contribution tends to satu- 
rate. Therefore, the T c should decrease as it is also ob- 
served in the experimenti 24 i 25 

In conclusion b y solving the nonlinear Eliash bcrg equa- 
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FIG. 4: Calculated superconducting transition tempera- 
ture T c (a) and quasiparticle effective mass m*(b) for 
Lai_ a! Pr 2 :Os4Sbi2 as a function of Pr concentration x. The 
straight curve is a guid e to the eye. The experimental data 
are taken from Ref. [lo- 
tions we find that in Lai_ a; Pr, ! ;Os4Sbi2 the dominant 
quadrupolar scattering due to the specific Th crystal 
group symmetry of the lattice is responsible for the ob- 
served increase of T c as a function of Pr concentration. 
Our analysis suggests that a combination of conventional 
electron-phonon interaction together with pairing medi- 
ated by quadrupolar excitons yields strong-coupling su- 
perconductivity in the PrOs4Sbi2 system. 

We are thankful to G. Zwicknagl, A. Yarcsko, S. Bur- 
din, and T. Westerkamp for helpful discussions. 
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